Separation of quadrupolar and magnetic contributions to spin— lattice relaxation in the 

case of a single isotope 
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We present a NMR pulse double-irradiation method which allows one to separate magnetic from 
quadrupolar contributions in the spin-lattice relaxation. The pulse sequence fully saturates one 
transition while another is observed. In the presence of a |Am[ = 2 quadrupolar contribution, the 
intensity of the observed line is altered compared to a standard spin-echo experiment. We calculated 
analytically this intensity change for spins / = 1,3/2,5/2, thus providing a quantitative analysis 
of the experimental results. Since the pulse sequence we used takes care of the absorbed radio- 
frequency power, no problems due to heating arise. The method is especially suited when only one 
NMR sensitive isotope is available. Different cross-checks were performed to prove the reliability 
of the obtained results. The applicability of this method is demonstrated by a study of the plane 
oxygen ^^O (/ — 5/2) in the high-temperature superconductor YBa2Cu408: the ^'^O spin-lattice 
relaxation rate consists of magnetic as well as quadrupolar contributions. 



PACS numbers: 



I. INTRODUCTION 



The work presented in this paper has been motivated 
by the experience in condensed matter nuclear mag- 
netic resonance (NMR) experiments that quite often both 
magnetic and quadrupolar time dependent interactions 
are present causing spin-lattice relaxation. The question 
arises whether it is possible to deduce, directly from the 
experiment, the admixture of the two different contribu- 
tions to the overall relaxation. 

The literature contains mainly calculations of multiex- 
ponential magnetization recovery laws for the case of ei- 
ther purely magnetic-jOr purely quadrupolar fluctuations, 
with Andrew et al.U being the first to treat the case 
of a static quadrupolar perturbed Zeeman Hamiltonian 
(spin / = 3/2|, 5/2). These calculations were extended to 
higher spinsoTM, and to the case of a si^atic quadrupolar 
Hamiltonian eI^eI. MacLaughlin et al.u treated the case 
of a static quadrupolar Hamiltonian ij] — 0) with mixed 
fluctuations in a kind of perturbation expansion, whereas 
Regally presented, for this case, an exact solution in the 
limit of time approacbing zero. 

In a previous studytil, we discussed the multiexponen- 
tial recovery for the case of a static quadrupolar per- 
turbed Zeeman Hamiltonian in the presence of both mag- 
netic and quadrupolar fluctuations under the assumption 
that the spin-exchange coupling can be omitted and the 
eigenf unctions of the static Hamiltonian can be approx- 
imated by Zeeman eigenf unctions. We found that, in a 
surprisingly large region of the parameter space spread by 
the probabilities for magnetic and quadrupolar induced 
transitions, it is almost impossible, within experimental 
errors, to separate magnetic and quadrupolar contribu- 
tions to the relaxation. Instead, the "dominant" con- 
tribution determines the time evolution of the recovery 
law, i.e. the system can approximately be described by 
a single time constant, T1 . However, it is questionable 
whether this approximation is meaningful in the presence 
of mixed relaxation or whether it is more appropriate to 



describe the system by the separate transition probabili- 
ties. 

If the nucleus under consideration has two magnetic 
isotopes as in the case of copper (^'^Cu and ^^Cu), the 
admixture can be estimated from the ratio of the relax- 
ation times, Ti. Nevertheless, the precision needed for a 
reliable interpretation of this ratio is commonly underes- 
timated. 

In this publication we will present a method, which 
enables the experimentalist to separate the different con- 
tributions of the spin-lattice relaxation especially in the 
case where the element under consideration has only one 
NMR sensitive isotope. The method involves a special 
initial condition of the spin system which we call dy- 
namic saturation and whichJiad already been mentioned 
briefly in our previous workLiJ. 

The paper is organized as follows. The next section will 
introduce the theoretical background essential to under- 
stand the method. In Sec. Ill we will give the results 



including a discussion, and Sec. IV will show experimen- 
tal results of ^^O NMR in YBa2Cu408 including a more 
technical discussion of the experiment. 



II. BASIC RELATIONS AND MASTER 
EQUATION 

For simplicity and the reader's convenience, we repeat 
part of our treatment presented in Rcf. Ell. The starting 
point is the following Hamiltonian: 



"^tot — ^0 



n^it), 



where Tio = 'Hz + ^q describes the time-independent 
(or "static") Hamiltonian which comprises the Zeeman 
interaction, TYz, with the external magnetic field and the 
quadrupolar interaction, Wq, with the internal electric 
field gradient (EFG) tensor. 7ii{t) takes into account 
fluctuations; it is the sum of a magnetic and a quadrupo- 
lar contribution: 



Hl(t) — Timngit) + Tiquadit), (1) 



where 



Wquad(t) = 47(1^ E Mt)T2k{I). 

Here, I is the nuclear spin operator, h{t) is a fluctuating 
magnetic field, Vk(t) is a component of the flu|etiiating 
EFG, and T2k{I) are spherical tensor operators EJ'EJ. 

In Eq. (|l|) , nuclear spin-exchange terms were omitted. 
If the quadrupolar splitting, due to TYq, is large compared 
to the nuclear spin-exchange coupling, the time evolu- 
tion of the spin-lattice relaxation proceeds by the direct 
coupling to the lattice. Cases where the nuclear sfj^ 
exchange terms are important are discussed in Refs. 

The relaxation of the spin system towards its thermo- 
dynamic equilibrium is described by the so-called master 
equation 



dt 



P{t)^W{P{t)-P{0)}. 



(2) 



Here, P{t) is the population vector of the different en- 
ergy levels with P(0) being the equilibrium value. The 
relaxation matris, W, is, in second order perturbation 
theory, given bylla 

Wo,f3 "^^ ^ J dT exp(ic^,^r) {a\ni (r) |/3) (/3|Hi (0) |a) 

— oo 

where \a), |/3) are eigenstates of Hq and ujaf3 — 
((ajTiolo^) — {l3\'Ho\/3))/fi are transition frequencies. En- 
semble averages are denoted by (...). 

As long as the eigenfunctions of Ho can be approxi- 
mated by the eigenfunctions of a Zeeman Hamiltonian, 
i.e. ||7iz|| ^ II^qII: the relevant relaxation matrix terms 
for magnetic and quadrupolar relaxation are given as fol- 
lows: 



W, 



ap 



Jicj^0)-{\{a\I+m' + \{a\r\f3}\'} 



+ \{a\ri,+IJ-\f3)\''} 
W/„7^'2^J(2)(,,„,).{|H(/+)2|^)|2+ 

+\{a\{i-rm'} . 

The J's are the spectral densities of the fluctuating fields: 
J {to) — —^ dr exp{iu}T) [ft,-|_,ft,_] 

^ J-oo 

j(i,2)(^) = f^\ p drexp(zc.T)[V+i,2,V-i,2] 



with [A, B] = (1/2) (A(t)B(O) + B{t)A{Q)), h± ^ K ± 
ihy, V±i = V^z ± iVyz, and V±2 = ^{V^x - Vyy) ± iV.^y. 

If TYz and Hq are of similar magnitude, the situation 
is more complicated. The case of purely magnetic fluc- 
tuations, for. ||?iz|| ~ II^qIIi has been treated by various 
authorsEiy. 

In this paper, we will deal with the case \\Hz\\ ^ II'^qII 
and we make the additional assumption that the spectral 
densities can be approximated by a single value. This 
means that the inverse of the correlation time, t^^, of 
the fluctuating fields is large compared to uja^, that is 
i-^apTc <C 1. One then obtains: 



J(^) 



J(0) =: W 
j(i.2)(o)^: ly^. 



and the resulting transition probabilities become 
W^J^^, ^W{I + m){I-m+l) 



W 



quad,l 



m—i-m — l 



w. 



w 



quad, 2 



m—^m — 2 



W2 



(2m - 1)2(7 - m + 1)(J + m) 

2/(27-1)2 
(7 + m) (7 + rn - 1) (7 - 771 + 1) (7 - m + 2) 
27(27-1)2 ■ 



Our calculations were performed in the high-temperature 
limit, i.e. hiOap <^ ksT, so that a further simplification 
takes place: Wa^p ~ Wp^a- Fig. |l| sketches the vari- 
ous transition probabilities which are possible for a 5/2 
spin system. We assume the spacings between the lev- 
els to be sufficiently unequal to suppress spin-exchange 
transitions. 

To solve the master equation, Eq. yj^, it is convenient 
to introduce some abbreviations. The population of level 
777 is Pm and we define the difference in population be- 
tween adjacent levels by Pra+1/2 = Pm+i — Pm', the equi- 
librium value of this difference is 779 — Pm+i(0) — Pm(0). 
The deviation of the population difference from its equi- 
librium value is denoted by Nm+1/2 = Pm+1/2 ^ ^0; the 
values iV„j+i/2 form the vector AT. 

Given the transition probabilities as shown in Fig. n, 
we can write down, in compact form, the following "re- 
duced" master equation for AT: 



— AT = RN, 

dt 



(3) 



where R is the reduced relaxation coefficient matrix. The 
solution of Eq. (H) is of the form 



^J-(^) = E [(E^) ' ^(0)] ^ E,j exp(tA,), 



(4) 



where A; and E are the eigenvalues and the eigenvector 
matrix of R, respectively, and E-^ denotes the transposed 
of E. N{0) is the vector describing the initial condition 
of the spin system into which it has been brought during 
a certain preparation period. 

Once the Nj{t) are known, the corresponding time de- 
pendent magnetization, M(t), is obtained: 



M(t) = Af (oo) 

and the a^ are given by 
1 r. 

Oi = 

no 



1 - y^Qiexp(tAi) 



E^ 



Ar(0) 



-Eli I 



(5) 



(6) 



where the index j refers to the hne which will be ob- 
served, e.g. the central transition. Usually, the irradiated 
line and the observed line are the same. 
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FIG. 1. Transitions between the spin energy levels effected 
by magnetic and quadrupolar spin-lattice relaxation pro- 
cesses for / = 5/2. 

We will now consider an experiment, which we call dy- 
namic saturation, with a special initial condition; this 
pulse sequence will allow us to disentangle quadrupo- 
lar from magnetic contributions in the spin-lattice re- 
laxation. We saturate a selected line, q, for instance by a 
long comb of pulses in such a way that the comb length, 
Ttot, is much larger than 1/ min(VF, Wi, W2) and that the 
pulse spacing, Tcd, within the comb satisfies the condi- 
tion 5T2 < Tcd < l/mSix{W,Wi,W2). This situation 
contrasts with an adiabatic manipulation of the spin sys- 
tem where, with the spin system initially in equilibrium, 
a short radio-frequency (RF) pulse is applied to one of 
the transitions. In the case of dynamic saturation, the 
initial condition must be calculated since the stimulating 
RF field causes transitions, with transition rate Prf, be- 
tween the levels q + 1/2 and g — 1/2. Thus, for calculating 
the initial condition vector, the rate equation ^) must be 
extended in the following way: 



dt 



(R+S)Ar + noP. 



s 



is a square matrix with all elements zero except 
Pii, S„„ — ~2Pi-f. P is a vector with all the ele- 



ments zero except Pq±i = P^t, Pq — —2P^{. For dynamic 
equilibrium, when dN /dt = 0, we have 



Ar(oo) 



-no(R+S) 



N{oo), which becomes the initial condition vector N{0) 
for solving Eq. (^), is calculated under the assumption 



that Prf ^ max ( W, VFi, 14/^2 )• The exact formulas for 
7V(0)/no are given in Appendix ^. 



III. TRANSITION ENHANCEMENT BY 
DYNAMIC SATURATION 

Let us deal, for the moment, with the special situation 
of pure magnetic relaxation, i.e. Wi,W2 — 0. In this 
case, after dynamic saturation, the Ar(0)'s for all spin 
values / > 1 take the form 

Ar(o)/7io = [o,...,o,-i,o,...,o] 

where —1 refers to the irradiated line. This equation re- 
flects the following behavior. On the time scale of the 
spin-lattice relaxation, Ti = 1/{2W), the new popula- 
tion differences are the same as in the case of thermody- 
namic equilibrium, except for the irradiated line. That 
means that a spectrum obtained by adiabatic manipula- 
tion {e.g. in a standard spin-echo experiment) is identical 
with the corresponding spectrum due to dynamic satu- 
ration, except for the irradiated line. This is not true 
anymore in the case of mixed or pure quadrupolar relax- 
ation as we will show now. 

The intensity of a specific transition which we observe 
will be denoted as follows. /ad|„i^,„_i is measured by 
an adiabatic pulse sequence as in the case of a standard 
7r/2 — TT spin echo experiment; /dyn|„Tlm_i I'^fers to the 
same transition, m -^ m — I, however, in the presence of 
dynamic saturation of the transition n — > n — 1. Given 
this notation, we define an enhancement factor by 



E2 



J in— >n— 1 



N{0) 
no 



m-1/2 



With the results of Appendix ^ we get, for instance, 
for a spin / = 5/2 system with the central transition 
being dynamically saturated and the inner satellite being 
observed: -^±3/2*_^±i/2 = 1 + Ms/Ci- The enhancement 
factor is one in the case of pure magnetic relaxation but 
it is nontrivial in the case of mixed or pure quadrupolar 
relaxation. 

That the enhancement factor is a nontrivial function, 
of the relaxation process was already noticed by PoundEZI 
who used it to show that ^^Na (/ = 3/2, 100% abun- 
dance) in NaNOa relaxes purely quadrupolar. 

Fig. y shows, for a spin / = 5/2 system with mixed 
relaxation, contour plots of the enhancement factor as a 
function of Wi/W and W2/W. Similar results, however 
less pronounced, are found for other combinations of 71 — > 
n — \ and m ^ m — 1. 

The characteristic results are as follows, (i) The en- 
hancement effect is less pronounced in the case of mixed 
relaxation as compared to pure quadrupolar relaxation. 
This makes it more difficult to detect the effect, although 
not impossible because of the very high time stability of 
modern spectrometers, (ii) There is always at least one 



transition with an appreciable enhancement factor [e.g. 
the cases (a) and (c) in Fig. ||], whereas the other tran- 
sitions are "depressed" [cases (b) and (d)]. This feature 
can be used for crosschecking the experiment if one is able 
to observe both transitions at the same time; this will be 
demonstrated in the next section, (iii) The enhancement 
factor depends only weakly on Wi/W since Wi connects, 
except for the (—1/2, l/2)-transition, the same levels as 
W does, which has no effect on the enhancement func- 
tion. Therefore, measuring the enhancement yields infor- 
mation only about the quadrupolar Am — 2 transitions. 
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FIG. 2. Contour plots of the enhancement for a spin 
I = 5/2 system in the case of mixed relaxation, (a) and 
(b) correspond to saturation of the central transition, (c) and 
(d) to saturation of the outer high-frequency satellite. 



IV. EXPERIMENTAL DETAILS 

We will discuss experimental details guided by 
our study of the high-temperature superconductor 
YBa2Cu4|Og; these investigations will be published 
elsewhereli3. 

The experiment was performed by using a combina- 
tion of two standard pulse spectrometers together with 
a magnetic field of Bq = 8.9945 T, {Bq \\ c). The reso- 
nant circuit was damped by a 12 17 resistor in order to 
achieve a broad frequency range. The resonance signals 
were obtained by a phase alternating add-subtract spin- 
echo technique similar to that one described in Ref. [l9| 



followed by Fourier transformation of the spin-echo. 

Each experiment consists of a certain combination of 
pulse sequences which are shown in Fig. g. To measure 
-^dynj we apply the saturation and the detection sequence. 
In the first sequence, dynamic saturation of the n —^ n — 1 
transition is achieved by applying pulses C and D. The 
spacing between these pulses. Ten, has to be much larger 
than T2 and much shorter than Ti, i.e. 5T2 < Tcd ^ 
Ti. The length of these pulses is chosen very large (20 
/is) in order to saturate the n — > n — 1 transition only. 
Furthermore, we change the phase between the C and D 
pulses by 90° to get rid of possible coherence effects. The 
total length of the saturation sequence, Tg, has to be of 
the order of the longest effective relaxation time. The 
detection sequence is the usual spin-echo tt/2 — tt pulse 
sequence. Here, we use very intense pulses (the tt/2 pulse 
length is about Ifis) to observe both the central transition 
and the high-frequency satellites in a single shot. 




Saturation-Sequence 



Heating-Sequence 



Detection- 
Sequence 

FIG. 3. Pulse sequence used in the experiment. Details are 
described in the text. 



In order to measure /ad, we apply, of course, the detec- 
tion sequence, but it must be supplemented by a heating 
sequence. Since pulsing heats the sample, a comparison 
of line intensities of different experiments (standard spin- 
echo and dynamic saturation) requires constant sample 
temperature. This is achieved by making the heating and 
the saturation sequence identical so that the total power 
absorbed by the sample is the same in either case. In 
other words, the total length of the pulse sequence, Tiot, 
is kept fixed (Ttot = 400 ms in our case). From our pri^i 
vious experiments on the isotope effect of the spin gapEHl 
we know that, at about 95 K, a constant temperature 
is achieved after running the heating sequence for about 
5 minutes. A combination of all three pulse sequences 
is used to cross-check our results; this will be discussed 
further below. 

To illustrate the method, we present ^^O spectra from 
YBa2Cu408 taken at T = 95 K, see Fig. |. This super- 
conductor contains, beside the apex oxygen, plane oxy- 
gen sites, 0(2) and 0(3), in the Cu02 plane, where super- 
conductivity takes place, and a chain oxygen site, 0(1 
Here, we are only interested in the plane sites. Fig. 
a) shows the "'^''O spectrum as obtained by Fourier trans- 
form of the standard echo. All central transitions coin- 
cide, all 0(1) satellites nearly coincide, while inner and 



outer 0(2,3) satellites are well separated. The splitting 
of these satellites is due to the orthorhombic symmetry 
(a = 3.8411 A, 6 = 3.8718 A) of the crystal. 
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FIG. 4. a) ^^O central transitions and high-frequency 
satelhtes of plane oxygen, 0(2) and 0(3), and chain oxy- 
gen sites, 0(1), in YBa2Cu40g. b) Spectrum obtained after 
dynamic saturation of the central transition. Bottom: Differ- 
ence of b) and a) spectrum. 

Fig. ^ b) presents the ^^O spectrum after dynamic sat- 
uration of the central transition and Fig. c) gives the 
difference between the "saturation" and the "standard" . 
Because of the short A and B pulses, all transitions can 
be observed. Obviously, the central transition is satu- 
rated, note the "negative" intensity. Whereas the cen- 
tral transition of the plane oxygen is totally saturated, 
the central transition as well as the satellites of chain oxy- 
gen are not. This is due to the fact that the chain oxygen 
nuclei have a considerably shorter Ti than the plane nu- 
clei and that the pulse sequence was optimized for plane 
oxygen. That this is true could be proven by the sym- 
metric experiment where we dynamically saturated the 
outer high-frequency satellite. 

The most important feature is the remaining positive 
intensity of the difference spectrum at the position of the 
inner high-frequency satellite. According to the discus- 
sion above and, in particular, the contour plot of Fig. 
(a), this intensity enhancement clearly shows that there 
is a quadrupolar contribution to the relaxation of the 
plane oxygen nuclei. This conclusion is supported by the 
fact that the intensity at the position of the outer high- 
frequency satellite is almost zero in agreement with the 
contour plot of Fig. 13 (b). 

The amount of quadrupolar admixture to the overall 
spin-lattice relaxation is determined as follows. Experi- 



mentally, the intensity enhancement of the inner satellite 
(when the central transition is dynamically saturated) 
is 1.13(2) which then results in a ratio W2/W = 1.4(3) 
according to Fig. || (a). If we saturate the outer high- 
frequency satellite, the intensity enhancement of the in- 
ner satellite is 1.04(2), leading to a ratio of W2/W = 
0.7(4) according to Fig. (c). The weighted average is: 
W2/W = 1.15(25). 

It is more difficult to determine the ratio Wi/W since 



the enhancement factor E 



1/2^-1/2 

±3/2-^±l/2 



depends only weakly 



on Wi. However, we can estimate this ratio indirectly, 
to be shown below; we found that Wi < W/3. 

To cross-check the results, we also performed an ex- 
periment with the so-called gradual saturation sequence 
which involves the application of all three pulse sequences 
as they are shown in Fig. ||. Tg is the duration of the sat- 
uration sequence. By I^ we denote the intensity of a line 
in case of gradual saturation while /q is the intensity of 
the line in the standard spin-echo sequence (including 
the heating sequence). The gradual intensity enhance- 
ment is then defined as Is/Iq- In Fig. g|, we have plotted 
this enhancement for the central transition (top figure) 
and for the outer satellite (bottom) as a function of T,. 
Bullets refer to inner satellites (top and bottom), open 
circles to outer satellites (top) , and triangles (bottom) to 
the central transition. 

For short Tg, the inner satellites (denoted by bullets) 
are strongly enhanced as expected, since in the case of an 
ideal 7r/2-pulse (adiabatic manipulation) at the central 
transition or outer satellite, respectively, Is/Iq should 
reach the value 1.5. Due to spin-lattice relaxation pro- 
cesses this value decays towards a limit which is (for Ts 
reaching a value corresponding to dynamical saturation) 
one in the case of pure magnetic relaxation and differ- 
ent from one in the case of mixed or pure quadrupolar 
relaxation. 

The response of the outer satellite and the central 
transition enhancement is retarded, since spin-lattice re- 
laxation processes need some time for "pumping" these 
transitions. This explains the enhancement maximum 
around Ts = 40 ms, before the enhancement starts to 
diminish towards one in the limit of dynamic satura- 
tion {Ts ~ 300 . . . 350 ms). In this limit, that is for 
Ts -^ 1/ mm{W,Wi,W2), we have, according to Sec. Ill 
and Fig. |, Is/Io ^ E. 
cussed above. 



1/2^-1/2 



This is the result dis- 
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FIG. 5. Top: Gradual saturation of central transition. Bul- 
lets (open symbols) represent the intensity enhancement of 
the inner (outer) high-frequency satellites. Bottom: Grad- 
ual saturation of outer high-frequency satellite. Bullets (tri- 
angles) represent the intensity enhancement of the inner 
high-frequency satellite (central transition). 

We are now able to estimate the ratio Wi/W. Is/Iq of 
the inner satellite decays very fast within the first 100 nis 
which corresponds to T^^ — 101,(5) ms at T = 95 K which 
we had measured previouslyEil by the standard inver- 
sion recovery method. However, dynamical saturation is 
achieved only at around 300. . . 350 ms. Therefore, a slow 
relaxation rate must be involved. Because W2/W « 1, 
we conclude that this slow relaxation rate has to be Wi, 
i.e. Wi < W/3. This clearly shows that the spin-lattice 



relaxation process contains a strong quadrupolar contri- 
bution which could not be detected otherwise so far. The 
discussion of the origin of this effect, which does not arise 
from phonons, as well as its temperature dependence will 
be given elsewheretS. 



V. SUMMARY 

We have presented a pulse double-irradiation method 
which allows one to separate magnetic from quadrupolar 
contributions in the spin-lattice relaxation. The pulse se- 
quence fully saturates one transition while another is ob- 
served. The clue is that the observed transition changes 
its intensity if and only if a | AttiJ = 2 quadrupolar contri- 
bution is present; the change is monitored with respect 
to a standard spin-echo experiment. We calculated an- 
alytically this intensity change for spins / = 1,3/2,5/2, 
thus providing a quantitative analysis of the experimen- 
tal results. Since the presented pulse sequence takes care 
of the absorbed radio-frequency power, no problems due 
to heating arise. The method is especially suited when 
only one NMR sensitive isotope is available. Different 
cross-checks were performed to prove the reliability of 
the obtained results. 

The applicability of the method is demonstrated by a 
study of the plane oxygen ^^O (/ — 5/2) in the high- 
temperature superconductor YBa2Cu408. We showed 
that the spin-lattice relaxation rate consists of magnetic 
as well as quadrupolar contributions. 
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APPENDIX A: ANALYTICAL FORMULAE FOR 
INITIAL CONDITION VECTOR 



N{0) _ M 

m C 



1. Spin 1=1 

^ = 2W2, C^'2W + Wi + 2W2 



2. Spin 1=3/2 

(a) Dynamic saturation of the central transition: 
N{0) 



no 



= [^ -1 ^1 



^i^W2, C = 3W + Wi + W2 



(b) Dynamic saturation of the satellite: 



-/V(O) ^ . Ml M2, 

"0 c c 



^1 ^ W2(3W + Wi + W2) 

M2 = -VK| 

C == 1214^2 + AWWi + IOWW2 



2W1W2 + wi 



+IO56WWIW2 + —W?W2 
5 



-27UWW,W^ 



392 



WiW^ 



5 ^ ^ 



16240H/2VFn^ 



98014^ VF,'' 



91^0" 



3. Spin 1=5/2 

(a) Dynamic saturation of the central transition: 



N(0) 
no 



rM4 M5 



M5 M4i 



Ci'Ci' 'Ci'C 



-9M^o2 



M4 — -^"2 
/i5 = 4514^2 (lOVK + 2Wi + W2) 
(1 = 4000T4^^ + lOOOVFH^i + 404^^ + IIOOVKVK2 
+ 160T4^iT4^2 + 45W|. 

(b) Dynamic saturation of the inner satellite: 

N(0) ^ rM6 _, Me M7 M8n 

no C2 C2 C2 C2 



^g = VK2 (800014^^ + 20001^^1^1 + 80T4^VFi^ + 

+380014^214^2 + 72014^14^1 W2 + 16T4^f 14^2 + 
+44014^14"! + 4614^1 14^1 + 914^|) 

M7 = -914^1(1014^ + 2VKi + 142)2 

Ms = 914^|(10VK + 214^1 + 14^2) 

C2 = 8000014^'' + 3600014^^14^1 + A^QQW^Wf + 

+16014^14^f + 820014^2|y2 ^ 4600014^^14^2 + 
+lQ%00W^WiW2 + \mQWWlW2 + 32^31^2 



-2O6OWW1W2 + lOSW^Wi + 5301F14f 



+641^114^1 



914^1 



(c) Dynamic saturation of the outer satellite: 
N{0) J _, M9 Mio Mil Mi2i 



no 



[-1, 



Cs ' Cs ' Cs ' Cs 



^9 = VK2(8000VF^ + 2000VK2^i + 8014^14^1^ _^ 

+3800W^W2 + 72C1WW1W2 + IQWIW2 
+UOWWi + AGWiWi + 9Wi) 



Mio = -14^1(80044^2 _^ 2004FVFi - 

+2204F4F2 + 32Wi4y2 ^ 
^11 = 9Wi{10W + 2Wi + W2) 

M12 = -914^1 
Cs = 128000VF^ + 3840041^^14^1 



- 814^1^ ^ 
914^1) 



■ 28804^2^2 _^ 
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